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The BBGKY hierarchy of equations for a particle interacting with an ideal gas is investigated. 
Principal properties of its solutions are disclosed, as exact identities which connect probability 
distribution of path of the particle, its derivatives in respect to gas density and irreducible many- 
particle correlations between gas molecules and the path. They show that all the correlations always 
give equally important contributions to evolution of the path distribution, and therefore the exact 
theory does not reduce to the classical kinetics even at arbitrary small gas density. 
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1. One of creators of the modern probability the- 
ory A. Kolmogorov underscored [l| that in it " the con- 
cept of independence of experiments fills most important 
place^^ and ^^correspondingly one of most important ob- 
jectives of philosophy of natural sciences " is " clearing- 
up and refinement of those prerequisites under which one 
can treat given phenomena as independenV\ Recall that 
in the probabihty theory some random phenomena or 
quantities A and B by definition are independent if 
their probability distributions are independent, that is 
P{A,B) — P{A) P{B) However, in natural sciences 
the independence of phenomena A and B is thought as 
absence of cause-and-cffect relations between them, that 
is an influence of one to another. Whether independence 
in this usual sense does mean independence in the sense 
of the probability theory? 

From viewpoints of common sense and philosophy, cer- 
tainly does not mean. Merely because A and B which 
do not directly influence one on another nevertheless both 
can be parts of a same another random event and thus 
turn out to be indirectly related. 

From the scientific point of view, it is natural to bring 
the same question to the statistical mechanics. One of 
creators of modern theory of dynamical systems and sta- 
tistical mechanics N. Krylov thoroughly analyzed it 
and also came to the negative answer: he concluded 
that opinions that phenomena which are "manifestly 
independent" should have independent probability distri- 
butions" , and the like, are nothing but prejudices" 0. 

Especially Krylov pointed [2] to firmness of such prej- 
udices 0]. Only it explains why the molecular chaos hy- 
pothesis put forward by Boltzmann many years ago \^ 
until now dominates kinetics although never was logically 
substantiated And why N. Bogolyubov, when he ob- 
tained ^ an exact hierarchy of evolution equations for 
s -particle distribution functions, straight away truncated 
his equations at s = 2 thus reducing it to the Boltzmann 
equation. 

Undoubtedly, molecules of sufficiently rarefied gas are 
independent in usual sense since almost surely have noth- 
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ing common in the past. Nevertheless they can be essen- 
tially dependent in the sense of the probability theory. 
This is quite understandable [7] (or you may see [1]). 
The absence of common causes of colliding particles in 
the past means, for each of them, absence of any back 
reaction of the gas to its past collisions. Therefore arbi- 
trary long fluctuations in relative frequency of collisions 
are allowable @. These fluctuations just play the role 
of aforesaid random events producing indirect statistical 
interdependencies between pairs (or groups) of molecules 
capable of being participators of one and the same colli- 
sion (or a cluster of successive collisions). 

As the consequence, P{A, B) ^ P{A) P{B) where 
P{A) is probability of flnding a molecule at (phase) point 
A and P{A, B) is probability of flnding simultaneously 
two molecules at points A and B . At that, relaxation of 
one-particle distribution P{A) is determined by pair cor- 
relation P[A, B) - P{A) P{B) . Relaxation of the latter 
just similarly always (regardless of the gas rarefaction) 
is determined by three-particle correlation. And so on 
up to infinity. Since during time interval t a molecular 
undergoes ^ t/r collisions (with r being characteristic 
free- flight time), a correct description of gas evolution 
over this interval requires taking into account s-particle 
correlations with at least s ^ t/r . Hence, in practice 
the whole hierarchy of equations deduced by Bogolyubov 
[1] is necessary. 

In work (and additionally or instead in Q) and 
in work approximate solutions to this hierarchy or, 
in other words, the Bogolyubov-Born-Green-Kirkwood- 
Yvon (BBGKY) equations were suggested for the 
problem about random wandering of a test molecule, and 
explanations were expounded why the Boltzmann's hy- 
pothesis is wrong. The aim of the present communication 
is to prove the statements of preceding paragraph with- 
out any approximations. At that we will strengthen the 
proof and besides simplify it due to replacing the usual 
gas by ideal gas whose molecules interact with the test 
molecule only but not with each other. 

2. We want to consider thermal random motion of 
a test molecule (TM) in thermodynamically equilibrium 
gas, at that specifying its position R(t) at some initial 
time moment i = : R(0) = Rq . 

Let P and M denote momentum and mass of TM, 
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m, Tj and Pj {j = 1,2,...) denote masses, coordinates 
and momenta of other molecules, $(r) is (short-range 
repulsive) potential of interaction between any of them 
and TM, and n is gas density (mean concentration of 
molecules) . At arbitrary time t > , full statistical 
description of this system is presented by the chain of 
{k + l)-particle distribution functions {k = 0,1,2,...): 
-fo(i, R; P| Ro ; n ) which is normalized (to unit) den- 
sity of probability distribution of TM's variables, and 
Ffe(i,R,rW,P,p«|Ro;n) (where rW = {ri...rfc}, 
pC^) — {pj... Pj, } ) which are probability densities of find- 
ing TM at point R with momentum P and simultane- 
ously finding out some k molecules at points Yj with 
momenta pj . A rigorous definition of such distribution 
functions (DF) was done in [gI]. In respect to the co- 
ordinates Tj they are not normalized, but instead (as 
in [^) obey the conditions of weakening of inter-particle 
correlations under spatial separation of particles. Sub- 
ject to the symmetry of DF in respect to Xj — {r^, Pj} 
these conditions can be compactly written as follows: 
Fu Fk-iG„,{pk) at Tfe ^ oo , where G™(p) is 
the Maxwell momentum distribution of a particle with 
mass TO . 

The enumerated DF satisfy a standard chain of the 
Bogolyubov equations [6]: 



dt 



[Hk,Fk]+n— $'(R-rfe+i)i^fe+i 
dP Jk+i 



{k = 0, 1, ... ) along with obvious initial conditions 



Fk 



fc t=o 



,5(R-Ro) exp(-iJfe/r) 



= (5(R - Ro) Gm(P) nLi E{v, - R) G,„(p,) , 



(1) 



(2) 



where is Hamiltonian of subsystem " k molecules 
+ TM", /^... = //... drfcdpfc, [...,...] means the 
Poisson brackets, <i>'(r) ~ V<l'(r) , and E{y) — 
exp[— $(r)/T] . Notice that TM can be considered as 
a molecule of non-uniformly distributed impurity, and 
equations ((J) are identical to the equations of two- 
component gas in the limit of infinitely rare impurity, 
when the main component is spatially homogeneous and 
thermodynamically equilibrium. 

Equations ([T]) together with ^ unambiguously deter- 
mine evolution of Fq and eventually probability distribu- 
tion of TM's displacement R — Rq . These equations will 
become more clear if we make a linear change of DF Fk 
by new functions Vk with the help of recurrent relations 
as follow: 

Fo(t,R,P|Ro;n) = K)(i, R, P| Ro; ") , 



Fi(i,R,ri,P,pi|Ro;n) = 
Fo(t,R,P|Ro;n)/(ri-R,pi)- 
+ Vi(i,R,ri,P,pi|Ro;n) , 



(3) 



i^2(i,R,r(2\P,p(^)|Ro;") = 
= Vo{t, R, P|Ro; n)/(ri - R, Pi)/(r2 - R, pa) + 
+ Vi{t, R, n, P, pi I Ro; n) /(ra - R, pa) + 
+ Vi{t, R, r2, P, P2I Ro; n) /(n - R, pi) + 
+ y2(i,R,r(2),P,p(2)|Ro;n) , 

and so on. 

Apparently, from viewpoint of the probability theory, 
Vk represent a kind of cumulants (semi-invariants), or cu- 
mulant functions (CF). It is important to notice that zero 
values of these CF would mean that all conditional DF 
of gas, Fk/Fo, are independent on initial position Rq 
of TM and thus on its displacement R — Ro . This fact 
makes visible an interesting specificity of CF Vk : they 
are irreducible correlations of not only current dynamic 
states of TM and k gas molecules but also correlations 
of all them with total previous TM's displacement. 

In terms of the CF the BBGKY hierarchy acquires a 
more complicated tridiagonal structure (we omit uninter- 
esting algebraic details): 

^ = [Hk,Vk] - r'^+O^fe+i + 



PfejG„(pfe)^'(rfe-R) 



d 



MT dP 



Vk-i (4) 



Here E'{r) — WE{r) , and Pkj symbolizes transposi- 
tion of the pairs of arguments Xj and Xk ■ On the other 
hand, initial conditions ^ and the above-mentioned con- 
ditions of weakening of correlations Q take very simple 
form: 



yo(0,R,P|Ro; n) =,5(R-Ro)Gm(P) , 
T4(0,R,rW,P,pW|Ro;n) = , 
14(i,R,rW,P,pW|Ro;n)^0 at r, ^ o 



(5) 



(1 < j < k ). Thus, as it should be with cumulants, CF 
Vk disappear under removal of even one of molecules. 

From these equations and initial conditions (as well 
as physical reasonings) it is clear that the reduction to 
zero in ([5]) realizes in an integrable way, so that integrals 
Vk = Ik+i ^k+i are finite. Let us consider them. By ap- 
plying the operation J^, to equations (|4]) one easy obtains 
equations 



-^^[Hk,Vk]+n-l^n^-rk,.)Vk.. 



d_ 
dP 



f $'(R-rfc+i)Ffe+i 

Jk+l 



(6) 



-T^ PkjGrniPk)E'{rk-R) 
i=i 



d 



MT dP 



V 



k-l 



where f{r,p)^E{r)Gmip), 



(with k = 0, 1, ...). Because of ([S]) initial conditions to 
these equations are zero: Vk{t = 0) = at any k . 
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Now, in addition to Vk , let us consider derivatives 
of CF in respect to the gas density, = dVk/dn . It 
is easy to see that differentiation of ^ in respect to n 
yields equations for the V^' which exactly coincide with 
^ after changing there Vk by . Besides, in view 
of ([5]), initial conditions to these equations again all are 
zero: V^{t = 0) = at any A: > . These observations 
strictly imply exact equalities = Vk , or 

1^ H(t,R,rW,P,pW|Ro;n) = (7) 
= / V^fc+i(t,R,r(^+i),P,p(^+i)|Ro;n) 

Jk+l 

This is main formal result of the present paper. 

3. The result ^ contains the proof promised in Sec.l. 
Indeed, equalities ([7]) show, firstly, that all the many- 
particle correlations between gas molecules and past dis- 
placement of test molecule (TM) really exist, i.e. dif- 
fer from zero. Secondly, all they have roughly one and 
the same order of magnitude. For instance, if comparing 
their integral values, due to ^ we can write, in natural 
dimensionless units, 

n''J...JJVkdP = n'=Fj"^(t,A;n) ^ CkVo{t,A;n) , 

where Voit,A;n)= / Vo(t, R, P|Ro; rfP is probabil- 
ity distribution of the TM's displacement A = R — Rq , 
V^'-'^'' (i. A; n) = d'^Voit, A;n)/dn'' are its derivatives in 
respect to gas density n , and Ck some numeric coeffi- 
cients. Hence, all the correlations are equally important, 
and none of them can be neglected if we aim at knowledge 
about true statistics of TM's random walk. 

For more details, let us suppose that (s -I- l)-particle 
correlation is so insignificant that one can assign Vs — 
in ([4]). At that, according to (H])-©, all higher-order 



correlations also will be rejected. Then, obviously, ac- 
cording to ([7]), distribution Vo(t, R, P|Ro; n) and thus 
Vo(i, A;n) must depend on n definitely as an (s — 1)- 
order polynomial. But, from the other hand, distribution 
Vq what follows from the truncated chain of equations 
^ certainly is absolutely non-polynomial function of n . 
With taking into account that equalities ([7]) do express 
exact properties of solutions to ([l])-® we see that very 
deep contradiction is on hand. 

This contradiction clearly prompts us that truncation 
of the BBGKY hierarchy leads to qualitative losses in its 
solution. 

Some possible losses already were characterized in Q 
and (and fir stly even much earlier in ^2*1) and in part 
filled up in 0, • Therefore here we confine ourselves 
(continuing 5-th paragraph of Sec.l) by remark that cut- 
ting of the (s -|- l)-particle correlation means cutting of 
s-th and higher statistical moments of fiuctuations in rel- 
ative frequency of TM's collisions with gas molecules (in 
other words, fluctuations in diffusivity of TM [l2|)- At 
s = 2 these fluctuations are completely ignored, and 
such truncated equations ^ yield a closed equation for 
Vb(i, R, P|Ro; which is equivalent to the Boltzmann- 
Lorentz equation [ll[. 

It is necessary to emphasize that above reasonings, as 
well as the exact relations are indifferent to a degree 
of the gas rarefaction. Consequently, one can state that 
the Boltzmann-Lorentz equation (moreover, all the clas- 
sical kinetics including the Boltzmann equation and its 
generalizations) does not represent a (low-density) limit 
of the exact statistical mechanical theory. The conven- 
tional kinetics is only (more or less adequate or carica- 
ture) probabilistic model of exact theory. Of course, in 
the latter also molecular chaos does prevail. But here it 
is much more rich, even if speaking about rarefied gas, 
and does not keep within naive probabilistic logics. 
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